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ABSTRACT 



Context. The Carina region is an excellent astrophysical laboratory for studying the feedback mechanisms of newly born, very 
massive stars within their natal giant molecular clouds (GMCs) at only 2.35 kpc distance. 

Aims. We use a clumpy PDR model to analyse the observed intensities of atomic carbon and CO and to derive the excitation 
conditions of the gas. 

Methods. The NANTEN2-4m submilhmeter telescope was used to map the [Cl] ^Pi -^ Po, '^^2 -^ Pi and GO 4-3, 7-6 hues in 
two 4' X 4' regions of Carina where molecular material interfaces with radiation from the massive star clusters. One region is 
the northern molecular cloud near the compact OB cluster Tr 14, and the second region is in the molecular cloud south of r;Car 
and Tr 16. These data were combined with ^''CO SEST spectra, HIRES/IRAS 60 ^m and 100 /im maps of the FIR continuum, 
and maps of 8 jj,in IRAC/Spitzer and MSX emission. 

Results. We used the HIRES far-infrared dust data to create a map of the FUV field heating the gas. The northern region shows 
an FUV field of a few 10^ in Draine units while the field of the southern region is about a factor 10 weaker. While the IRAC 
8 fj,in emission lights up at the edges of the molecular clouds, CO and also [Cl] appear to trace the H2 gas column density. The 
northern region shows a complex velocity and spatial structure, while the southern region shows an edge-on PDR with a single 
Gaussian velocity component. We constructed models consisting of an ensemble of small spherically symmetric PDR clumps 
within the 38" beam (0.43 pc), which follow canonical power-law mass and mass-size distributions. We find that an average 
local clump density of 2 10^ cm^^ is needed to reproduce the observed line emission at two selected interface positions. 
Conclusions. Stationary, clumpy PDR models reproduce the observed cooling lines of atomic carbon and CO at two positions in 
the Carina Nebula. 

Key words. ISM: clouds, ISM: structure, ISM: individual objects: Carina 



1. Introduction 

1.1. The importance of atomic carbon 

The fine structure lines of atomic carbon are strong 
coolants of the interstellar medi um (ISM) through out the 
universe: in Galactic clo uds (e.g. Jakob et al.ll2007l) . in the 
Milky Way as a whole ( Fixsen et al.lll999l ). as well as in 
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external galaxies ( Bavet et al. 20061: iKramer et al.ll2005[) 
up to high redshifts (jWeiss et al.ll2003h . Large-scale obser- 
vations of Galactic clouds show [G i] emission coincident 
with that of GO. This hol ds for giant molecular clouds 
(GMGs) (ISakai et al.l 2006), as wel l as fo r rather quies- 
cent dark clouds (jTatematsu et al.l 119991 ). It has there- 
fore been suggested that [Gi] is a more reliable tracer 
of t he molecular cloud inasses in external galaxies than 
GO (jPapadopoulos et al.ll2004[ l. However, all steady-state 
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models of photon-dominated regions (PDRs) predict that 
atomic carbon is formed in a surface layer by recombi- 
nation of 0+ and dissociation of CO. This appears to 
be in contradiction with atomic carbon emission which 
is observed even at large distances from UV sources. The 
dynamics of the ISM may play a decisive role, allowing 
to explain the di stribution of [Ci] with time depend ent 
chemical models (|Oka et al.ll2004HSt6rzer et al.lll997t ). 

However, stationary, but clumpy PDR models have 
been successful in explaining e.g. extende d [Ci] ^Pi —^ 
Pn and 1^00 2 - 1 e i nissio n in S 140 (jSpaansI llQQfil : 
Spaans fc van Dishoecla Il997t ). Inside a clumpy PDR, 
the UV radiation can penetrate much d eeper than i t 
would be possible ii i a ho mogeneous layer (|Boisselll990l ). 
Meixner fc Tielend ( 19931 ) showed that the line intensi- 
ties of all lines are enhanced compare d to a homogeneous 
model with the same average density. Ijuvela et al.l (|200l[ ) 
presented a model to apply Monte Carlo radiative transfer 
calculations to the results of magnetohydrodynamic cal- 
culations to account for a clumpy and turbulent structure. 
They also found that local density condensations lead to 
an increased cooling efficiency and reduced photon trap- 
ping. Another approach is to describe both the turbulent 
velocity field and the density s tructu re by their statisti- 
cal properties. iHegmann et al.l (J2007I ) calculate CO cool- 
ing rates based on a stochastic radiative transfer model, 
which accounts for density and velocity fluctuations with 
a finite correlation length. 

Here, we study the distribution of atomic carbon 
and warm CO in the Carina star forming region and 
test clumpy models using the KOS MA-r PDR model 
( Storzer et anil996HR611ig et al.ll2007[) to derive the exci- 
tation conditions of the gas. Spaans et al. model the inho- 
mogeneous structure with a Monte-Carlo code, assuming a 
fixed clump/interclump density ratio of 10. In contrast, we 
assume that all emission analyzed stems from the clumps. 
The influence of the interclump medium or a ha lo lead- 
ing to pre- shielding of CO has been discussed bv iBensch 
( 2006[) and lBenschet al.l (|2003l ) . who also used KOSMA-t, 
but is not considered here. To describe the clumpy struc- 
ture, we assume that the clumps are distributed follow- 
ing the canonical clump mass and mass-size distributions 
which have been found in molecular clouds. 



1.2. The Carina nebula 

The Carina nebula is a very prominent souther n star- 
formin g region at a distance of 2.35 kpc ( SmithI 
2006 a). This giant mo lecular cloud was first observed by 
Grabelskv et all (|l988l ) as part of the Columbia CO 1-0 
survey of the Milky Way. They derived a total mass of 
6.7 10^ Mq and an effective radius of 66 pc. In total 65 O- 
type stars, incl uding 6 of the 11 03-type stars known in 
the Milky Way (|SmitlTll2006"ah . produce a strong UV field 
and strong stellar winds which interact vigorously with the 
material in the surrounding molecular clouds. This high 



concentration of the earliest-type stars is unique in the 
Galaxy and may serve as a template for more extreme 
but more distant regions containing stellar super clus- 
ters like the central regions of NGC 253 or the Antennae 
(JBavet et al.ll2006HSchulz et al.ll2007l ). The Carina nebula 
extends more than 4 square degrees on the sky. It exhibits 
a peculiar morphological structure on all scales. Massive 
stars are being born in several regions within molecu- 
lar condensations, located mainly to the south-east and 
north-west of the cluster Trumpler 16, which is centered 
at TyCar, a luminous blue variable star. nCar is embed- 
ded in the Homunculus nebula (e.g. ISmithll2006a[ l and its 
luminosity has changed drastically and repeatedly in the 
past centuries, leading to strong variations of the FUV 
field impinging on the surrounding clouds. Its current lu- 
minosity of 5 10® LfT) inakes it one of the brightest celestial 



infrared objects (|Coxlll995l ). In the northern part of the 
Carina nebula, the stell ar population i s dom inated by the 
cluster Tr 14. Recentlv. lOberst et al.l (2006!) reported the 
detection of the [N ii] line at 205 /im at the radio peak of 
the Carina II H ii region using SPI FI at AST /RO. They 
comp ared their data with [C ii] data ( Mizutani et al.ll2004 . 
20021 ) and concluded that about a third of the [C ii] emis- 
sion stems from a low-density ionized medium. For a re- 
cent o verview of the Carina region, see ISmith fc Brookd 
()2007l) . 



Carina was mapped at ~ 3 ^ resolution in the C O 4-3 
and [C i] ^Pi -^ Pq lines by IZhang et all (|200lh using 
the AST/RO 1.7 m telescope. These AST/RO data show 
large-scale [Ci] emission, coextensive with CO 4-3. 

Here, we have used the 4 m NANTEN2 telescope to 
map ^^CO in the rotational transitions J =4-3, 7-6, and 
[Ci] in the fine structure transitions ^Pi—^Pq, and 
^P2 —^ Pi (henceforth 1-0 and 2-1) at sub-arcminute res- 
olutions in two 4' X 4' regions of Carina (Fig.[T]). 

One region to the north of ryCar covers the cloud inter- 
face with Tr 14. The second region covers another cloud 
interface to the south of r^Car and Tr 16. 

The northern molecular cloud to the west of Tr 14 
(Fig[T]) contains several sites of massive star forma- 
tion. Diffuse PAH 3.29 /J,m emission has been found by 
iRathborne et al.l (2002) to trace the sharp edge of the 
cloud. Peak intensities at the edge towards the Car I H ii 
region arc almost 10 times higher than in the other re- 
gions of Carina. Car I-E is interacting with th e front face of 



the CMC and creating a P DR seen edge-on (jBrooks et al 



20031 : iBrooks et all l200lL Carl-W manifests a second 



PDR seen face-on (cf. Fig. [2]). 



The CO 2-1 map of [c3 (|l995l ) shows a sharp 
rim of gas delineating the southern molecular cloud, lo- 
cated south of r/Car and the Tr 16 cluster. Several IRAS 
sources are located along the edge. A prominent peak in 
PAH 3.29 /im emis sion is the ultra comp act Hii region 
IRAS 10430 - 5931 ()Rathborne et al.ll2002h . 
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Fig. 1. Overview of the Carina clouds near jyCar lying at 0/0 10:45:03.59 -59:41:04.3 (J2000) at the center of the Tr 16 
cluster. Left: Estimate of the FUV flux x in Draine units at a resolution of ~ 1' and derived from HIRES/IRAS. 
Dashed contours range from 400 to 800 in steps of 200. Solid contours range from 10^ to 5 10'^ in steps of 10^. Right: 
MSX emission at 8 iim at a resolution of ~ 20". Contours are 1 to 10 in steps of 1 10~^ Wm~^sr~^ (cf. Figs. 7, 11 in 
iRathborne et al.ll2002l ). Note the diffraction artefacts near TyCar. The center position of the OB cluster Tr 14 is marked 
by a triangle. The two 4' x 4' regions mapped with NANTEN2 are delineated by boxes. Crosses in the left panel mark 
two cuts discussed in the text. 



Table 1. Physical parameters derived from the observed integrated intensities 7 in K kms ^ and intensity ratios Rcx)—I{CO 7- 
6)//(C0 4-3), i?c =-''([Ci]2-l)/7([Cl]l-0), i?cco=-^([Cill-0)//(CO 4-3) from Gaussian fits to the line profiles at two interface 
positions. The calibration error of the integrated intensities is estimated to be 15%. 
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2. NANTEN2 observations 

We used the new NANTEN2 4m telescope situated at 
4865 m altitude at Pampa La Bola in northern Chile 
to map two regions in Carina in CO 4-3 and 7-6, and 
in the two fine structure transitions of atomic carbon. 
Observations were conducted in 2006 between September 
21 and October 6 with a dual-channel 460/810 GHz re- 
ceiver installed for verifying the telescope submillimeter 
performance. Double-sideband (DSB) receiver tempera- 
tures were ~ 250 K in the lower channel and ^ 750 K 
in the upper one. The intermediate frequencies (IF) are 
4 GHz and 1.5 GHz, respectively. The latter IF allows si- 
multaneous observations of the CO 7-6 line in the lower 
and the [Ci] 2-1 line in the upper sideband. These two 
lines are observed simultaneously with one of the lines in 
the 460 GHz channel. As backends, we used two acousto- 
optical spectrometers (AOS) with a bandwidth of 1 GHz 
and a channel resolution of 0.37kms~^ at 460 GHz and 
0.21kms-i at 806 GHz. 



The pointing was checked regularly on Jupiter, 
IRC+10216, and IRc2 in Orion A. The applied correc- 
tions were always smaller than 20", and usually less than 
10". The atmospheric transmission was derived by mea- 
suring the atmospheric emission at the reference position. 
Carina was observed between ~ 32° and ~ 54° elevation. 
Spectra of the two frequency bands were calibrated sepa- 
rately, and sideband im balances were corr ected using the 
atmospheric model atm ( Pardo et al.ll2001l ). We used a ref- 
erence position at 10l'5a"09!3 -59?18:33"7 (J2000), - 45' 
to the north-e ast of ??Car, whi c h was observed to be free 
of emission bv lYonekura et al.l (20051). Observations were 
conducted on-the-fly (OTF), scanning in right ascension 
at a speed of 2.5"/sec and with a sampling interval of 
10". [C i] emission of the northern field was also observed 
scanning in declination. The reference position was ob- 
served at the beginning of each OTF scanning line, i.e. 
every 1.6 minutes in time. To complete a single scan of the 
4' X 4' map area took ~ 1 hour of total observing time. ON- 
sourcc integration times per position on a 10" grid range 
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between 24 sec (for CO 4-3 & 7-6 in the northern field) 
and 4 sec. The total observing time was 10 hours. 

The half power beamwidths (HPBWs) and main beam 
efficiencies (-Bcft) were determined fr om continuum cross 
scans on Jupiter ( Simon et al.ll2007l l. For the latter, we 
linearly i nterpolated the Jov ian brightness temperatures 
listed bv I Griffin et al.l (|l986h to the observed frequencies. 
The HPBWs deconvolved from the observed full widths 
at half maximum (FWHM), are 38'.'0 and 26'.'5, in the 
lower and upper receiver band, respectively. Beam effi- 
ciencies are 50% and 45%, respectively. The forward effi- 
ciency (fcff) was derived from sky dips and was found to 
be 86% for both frequency bands. The raw data taken at 
the telescope was recalibrated to the T^ antenna temper- 
ature scale and then multiplied by Fcs/BcS to scale to 
main beam temperatures. All data presented here are on 
the Tnih scale. We fitted and subtracted 2nd order poly- 
nomials from all spectra. In a few cases, we subtracted 
higher-order polynomials. 

During the test campaign in the second half of 2006, we 
were able to roughly meas ure the error beam amplitudes of 
the NANTEN2 antenna fiSimon et al.ll2007t ). Continuum 



measurements of the solar and lunar edges allowed to 
measure the extent and strength of a first compact er- 
ror beain_St.810 GHz following the method described e.g. 
bv iGreve et "ail (|l998h . We find a FWHM of 240" corre- 
sponding roughly to the average panel size of about 0.6 m. 
Only about 10-15% of the power is detected within this 
error beam which we attribute to panel misalignment. The 
power detected within a beam of ~ 30' was derived from 
lunar continuum scans using a lunar Rayleig h-Jeans tem- 



perature of 350 K fo r both frequency bands (jPardo et al 
2005HMangi!^ll993[) . Moon efficiencies are ~ 70% at both 



frequencies. In total, about 36% of the total power is de- 
tected in at least two error beams at both frequencies. 

As the spatial source structure in each velocity channel 
is convolved by the antenna diagram during observations, 
the error beams lead to pickup of extended emission in 
Carina. Main beam temperatures presented here are somc- 
what too high. Lin e profiles may be slightly distorted (cf. 
Bensch et al.ll2001l ). However, as the error beams are not 
yet completely characterized and may moreover be some- 
what time variable, we ignore their influence in the fol- 
lowing. To first order, [Ci] and CO emission stems from 
the same region. Line ratios will therefore be much less 
affected by the error beam compared to absolute intensi- 
ties. 



3. Complementary data 

3.1. ^^ CO 2-1 data 

Brooks et al. obtained complementary maps of ^^CO 2- 
1 at the Swedish-ESO SubmiUimetre Telescope (SEST). 
The telescope main beam efficiency at 220 GHz was 50% 
and its HPBW 25 ". The northern 4' x 5' SEST map 
(JBrooks et al.ll2003l ) covers the northern NANTEN2 field. 
The southern map (Brooks, Schneider, priv. comm.) cov- 



ers the southern NANTEN2 field. We convolved these data 
with a Gaussian kernel to a resolution of 38" to allow a 
comparison with the NANTEN2 spectra taken along two 
cuts through the interface regions. 

3.2. Dust and PAHs in Carina 

Figurc[T] gives an overview of a 25' x 25' region surround- 
ing rjCai, showing the Keyhole region in the center, and 
the northern and southern clouds. The OB stars of Tr 14 
and Tr 16 emit far-UV photons which are partly absorbed 
by the dust grains of the surrounding molecular clouds. 
The heated dust cools in the far-infrared via dust contin- 
uum emission. Gas heating is very inefficient; only a few 
percent at most of the absorbed energy is transferred to 
the gas (IHollenbach fc Tielens 1999 ). which cools via far- 
infrared and submillimeter line emission. Since the FUV 
photons play a decisive role in driving the cloud chemistry, 
determination of the FUV field is important. 

To derive a map of the estimated FUV fluxes (Fig.[T]). 
we obtained HIRES/IRAS 60 and 100 /im images at -- 1' 
resolution from the IPAC data center. These images at 
enhanced resol ution were created us ing a maximum corre- 
lation method ( Aumann et al.lll990h . The fllter properties 
of IRAS allow us to combine these two data sets to cre- 
ate a map of fa r-infrared intensities Jftr between 42.5 ur a 
and 122. 5^m (JHelou et al.lll988t [Nakagawa et aLlll998l ). 
Assuming that the FUV energy absorbed by the grains is 
reradiated in the far-infrared, we then estimate the FUV 
flux X (6eV < hv <13.6cV) impinging onto the cloud 
surfaces from the emergent FIR intensities: 



X/Xo = 47r /fir/ (ergs 



^cm ^sr ^) 



with Y in units of Xo 



19781: iDraine fc Bertoldi 



2.710- 
1996h . 



ergs cm 
Here, we assume 



^ ( Draind 
that 



heat ing by photons with hv < 6e V contributes a factor of 
~ 2 ( Tielens fc Hollenbachlll985r i and that the bolometric 
dust continuum inten sity is a factor of ~ 2 larger than 



/pTB (|Dale et al.l 120011 ). In this case, the two corrections 
cancel out. In case the clouds do not fill the beam, the de- 
rived estimate of the FUV field will only be a lower limit. 
Obviously, this method fails and the FUV field cannot be 
derived at positions without dusty clouds. 

The northern 4' x 4' region mapped with NANTEN2 
shows FUV field strengths varying between 2 10'^ in the 
north-east and 5 lO'^ in the south-west. The FUV field in 
the southern region is a factor 10 weaker, varying between 
4 1 0^ in the outs k irts an d 8 10^ near the center. 

Brooks et all (|2003l ) estimate the FUV field at the 
position of peak [C ii] emission in the northern cloud 
(Figs. 1113) 1 from the spectral types of the 13 brightest 
stars dominating the stellar content of Tr 14. They derive 
Go = 1.4 10^ in units of the Habing field, i.e. x = 8.2 10^, 
a factor 1.6 larger than the peak flux in the northern re- 
gio n, as derived f rom H IRES. However, the value derived 
bv iBrooks et al.l (|2003[ ) is a strict upper limit as it is as- 
sumed that the stars and cloud lie at the same distance 
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to the observer and attenuation by dust can be neglected. 
Moreover, the FUV field derived from the stars would be 
lowered by smearing to the resolution of the HIRES data. 
Figurc[l] shows contours of the FUV field together with 
an 8 iim MSX image. The MSX im age shows the emission 
from very small grains and PAHs (|Smithll2006a[ ) excited 
by the FUV field. Evidence that the PAH emission is stim- 
ulated by FUV absorption is reflected by the good corre- 
lation between the FIR continuum and the 8 /im emission. 




Fig. 2. Northern field: Contours of integrated CO 4-3 in- 
tensities (cf. Fig.[3|) overlayed on 8^m IRAC/Spitzer emis- 
sion in units of MJy/sr. C ar I-E and Ca r I-W a re two Hil 
regions described e.g. by iBrooks et al.l |2003[) (cf. Fig. [3] 
for more details). 




Fig. 3. Northern 4' x 4' field: Velocity-integrated maps of 
CO 4-3, 7-6, [Ci] 1-0, and 2-1 at a common angular res- 
olution of 38" (0.43 pc), integrated over the full velocity 
range of emission between —32 and — 5kms~^. Contours 
range between 10 and 90% in steps of 10% of the peak in- 
tensities, which are 398Kkms~^ for CO 4-3, 153Kkms~^ 
for CO 7-6, 59Kkms-i for [Ci] 1-0, and 48Kkms-i for 
[Ci] 2-1. The center position of the compact OB cluster 
Tr 14 is marked by a filled triangle outside the CO 4-3 
map. A dashed line and small crosses space d by 40" mark 



a cut through the peak of [Cil] emission (jBrooks et al 



120031 ) (large cross) and Tr 14. The circle marks the reso- 
lution and the interface position analyzed in Table[TJ 



4. Results 

4.1. Northern region 

4.1.1. Maps of integrated intensities 

Figure[2] shows an overlay of integrated CO 4-3 intensities 
on an 8 /xm image taken with the Infrared Array Camera 
(IRAC) on the Spitzer space telescopco- The 8/xm data 
reveal the complicated spatial structure of the region, 
consisting of several filaments and shell-like structures. 
The bright region in the south-west is the ionisation front 
Car I-W. Approximately ~ 1' to the east lies Car I-E (cf. 
Rathborne et al.ll2002HBrooks et alll200l[) . 

Maps of integrated CO 4-3, 7-6, [Ci] 1-0, and 2- 
1 emission (Fig. [3]) show the boundary of the north- 
ern molecular cloud (orientated north-south) towards the 
Trl4 cluster and its Hii region. For all four tracers, 
the integrated intensities are co-extensive and peak in 
the south-west, near —1 3V + 6.5^ i.e . near the peak of 
[Cll] emission studied bv lBrooks et al.l (|2003l ). The south- 
eastern interface region exhibits increased CO 7-6 and 



^ These post-BCD data were retrieved from the Spitzer 
archive, reproject ed, and merged to create Figs. [5] and [S] (cf. 



ISmithet aLll2005l l. 



[C i] 2-1 line intensities, indicating an increased amount of 
warm, dense gas. The morphology of the cloud resembles 
a ring centered on a region with less strong emission near 
— 12'/+6.'5. The northern half of the cloud shows weaker 
and more diffuse emission. 

The northern region exhibits a complex velocity struc- 
ture between —32 and —5 kms~^, which can be seen in the 
CO 4-3 velocity channel maps (Fig.|4]). Note that weak CO 
4-3 emission is also detected east of the north-south inter- 
face, near — 10'/-|-7.'5, at velocities between about —14.5 
and — lOkms^^. The CO 4-3 emission also corresponds 
weU with the ^^CO 2-1 S EST channel i naps at 25" reso- 
lution shown in Fig. 7 of iBrooks et al.l ( 2003[ ). The good 
corrcspondance shows that the velocity structure is largely 
due to individual clumps and not self-absorption of opti- 
cally thick 12C0. 

4.1.2. Spectra along a cut through the region 

The rich kinematical structure is further revealed by the 
spectra along a cut, shown in Fig.[5l from the [Cii] emis- 
sion peak through the cloud interface, towards Tr 14. This 
cut is ori entated the same as the cut presented in Figs. 1 
and 11 of IBrooks etlll (|2003f) . The spectra of CO 4-3, 7- 
6, [Ci] 1-0, 2-1, and "CO 2-1 (IBrooks et al.ll2003l) show 
several velocity components with large velocity gradients. 
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Fig. 4. Velocity structure of the northern field: CO 4- 
3 velocity channels of 1.5 kms^^ width. Contours range 
between 5 and 45Kkms^^ in steps of 5Kknis~^. 



~^,^:/!!;ff?;?'li~,^i^;i7!i{f?;Iili~'^;^!^f?;?'li,,J,7J:?;p!/?'*,,i 




Fig. 5. Spectra along a cut through the northern field con- 
necting Tr 14 and the [C ii] peak on a 40" grid. All data are 
at a common resolution of 38" and on the main beam scale. 
The velocity range is —40 to 5kms~^. The dashed vertical 
lines and the arrow mark the position of the — 24kms~^ 
velocity component presented in Table[T] 



At the interface, all tracers show a strong component at 
— 24kms~^. Further west, this component shifts to lower 
velocities and splits up into at least 3 components. The 
line centroids of the different tracers match very well. The 
peak temperatures are 40 K in CO 4-3 and 10 K in [Ci] 
1-0. Line intensities of CO and [C i] stay constant within a 
factor of 2 at the positions within the cloud; they increase 
slightly at the interface, before dropping rapidly towards 
the east. A weak component at —15 kms^^ seen in CO and 
^■^CO peaks to the east of the interface at (— 9f8,7f7). 

To give one typical example, Table[T] shows the in- 
tegrated line intensities and derived quantities of the 
— 24kms~^ component at the interface position, after 
Gaussian fitting. Tex is derived from the ^^CO 4-3 peak 
line temperature assuming optically thick emission and a 
beam-filling factor of 1. Total CO and C column densi- 
ties, given in lO^^cm"^, were derived from ^^CO 2-1 and 



[C i] 1-0 assuming optically thin emission at Tox, LTE, 
and a ^^CO vs. ^"^CO abundance ratio of 65. Total H2 
column densities, given in lO^^cm"^, and masses were 
derived f rom NjGO) assumin g an abundance ratio of 
8.510-5 (JFrerking et al.lll982l ). C/CO is the abundance 
ratio 7V(C)/7V(CO). Errors given for the column densities 
were derived by varying T^x and the integrated intensi- 
ties by ±15%. While the ^^CO 4-3 and 7-6 lines show 
FWHMs of 5.4 and 4.7kms-\ the [Ci] 1-0 and 2-1 line 
widths are smaller, 3.9 and 3.1 kms^^. ^^CO 2-1 shows a 
still narrower line width of 2.9 kms"^. We attribute this 
change of line widths to optical depth effects. The fitted 
line center velocities of all tracers agree within 0.5 kms"^. 
Line ratios of CO 7-6/4-3 and [Ci] 2-1/1-0 are 0.4 and 
0.7, which are typical values for the entire cut. 

The peak CO 4-3 temperature of 40 K translates, 
through the Rayleigh- Jeans correction, into a lower limit 
of the gas kinetic temperature of the emission zone of 
50 K. Optical depth effects or a beam-filling factor of 
less than 1 w o uld im ply higher gas kinetic temperatures. 
Brooks et al.l (|2003l) derived a dust temperature of 50 K 
from the 60 vs. 100 /xm IRAS flux ratio at the [C 11] posi- 
tion. 

The [Ci] 2-1/1-0 line ratio is a sensitive function of 
the [C i] excitation temperature. In the optically thin limit 
and assuming LTE, T^^ = 38.3K/ln[2.11/i?^\]. At the in- 
terface position, we observed a ratio of 0.72 ± 0.15, as- 
suming a calibration error of the ratio of about 20%. This 
translates into rex([Ci]) of only 35K(43K). 

The ^^CO 2-1 spectra (iBrooks et al.l l2003l) peak at 
15 K. At the interface, the [Ci] I-O/^^cq 2-1 ratio is 0.8. 

To derive a first estimate of the total CO column den- 
sities at the two interface positions (Table[T|), we used the 
^3 CO 2-1 integrated intensities assuming optically thin 
emission, LTE, and a ^^CO vs. ^"^CO abundance ratio of 
65. The abundance ratio is in accordance with the av- 
er age local ISM i2q/1 3q ^.^tio of 68 ± 15 recently found 
bv iMilam et al.l ( 2005[ ). which is valid also for Carina at 
about solar galacto centric distance. The inferred total op- 
tical extinction at the northern interface is 28 ±7 mag, us- 
i ng the canonical JV (Il2)/Av ratio of 9.36 10^" cm"^ mag 



(|Bohhn et al.lll978[ ). The C/CO abundance ratio is 0.21. 



4.2. Southern region 

4.2.1. Maps of integrated intensities 

FigurejBj shows an overlay of integrated CO 4-3 intensities 
on an 8 fim IRAC/Spitzer image of the southern NANTEN 
field. The 8/xm emission traces a sharp interface running 
from north-east to south-west. It brightens at the cloud 
edge, showing several elongated filaments and knots along 
the interface. It also shows a bright emission knob in the 
north-west, which probably corresponds to the IRAS point 
source 10430-5931. The 8 /im, emission here is weaker than 
in the northern field by about a factor 5. The IRAC emis- 
sion drops rapidly towards the bulk of the cloud in the 
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Fig. 6. Southern field: Contours of integrated CO 4-3 in- 
tensities (of. Fig.[7|) overlayed on 8 ^m IRAC/Spitzer emis- 
sion in units of MJy/sr. 




Fig. 7. Southern 4' x 4' field: Velocity-integrated maps 
of CO 4-3, 7-6, [Ci] 1-0, and 2-1 at a common angular 
resolution of 38", integrated over the full velocity range 
of emission between —20 and — 32kms~^. Contours range 
between 10 and 90% in steps of 10% of the peak intensities 
which are 114Kkms-i for CO 4-3, 35Kkms-i for CO 
7-6, 37Kkms-i for [Ci] 1-0, and 22Kkms-i for [Ci] 2- 
1. A cross X marks the position of IRAS 10430-5931 at 
ltfM5"^02!4 -59?47fl0'.'0 (J2000). Small crosses + mark 
positions along a cut through the interface region. The 
circle marks the resolution and the position analyzed in 
TableB 



south-ea st, similar to PAH emiss ion traced by the 3.21 fim 



image of iRathborne et al.l |2002l ) 



In contrast to the mid-infrared emission, CO 4-3 and 
[C i] integrated intensities (Fig.[71) peak south of the inter- 
face deep inside the cloud near (1', —8'). The CO 7-6 hue 
peaks near IRAS 10430-5931, indicating elevated temper- 




Fig. 8. Spectra along a north-south orientated cut 
through the southern field at Aa =0.'8 on a 40" grid. 
All data are at a common resolution of 38". The velocity 
range is —43 to — 13kms~^. Vertical lines mark a velocity 
of — 28kms~^. The arrow marks the position analyzed in 
TableE 



atures and densities. In contrast, the IRAS source is not 
prominent in [C i] . 

4.2.2. Spectra along a cut through the region 

Spectra taken along a north-south cut through the south- 
ern field show one component near — 27kms~"'^ with little 
variation in velocity (Fig.[8|). Near the interface. ^^CO line 
profiles are flat-topped indicating high optical depths and 
self-absorption. Table[l] lists the results of Gaussian fits to 
the spectra at the interface position (0f8,— 7f3). The fit- 
ted FWHMs are 4.3kms-i (CO 4-3), 3.9kms-i (CO 7- 
6), 3.4kms-i (i^CO 2-1), 2.1 kms'^ ([Ci]l-0), 2.4kms-i 
([C i]2-l). ^^CO line widths are a factor ~ 2 broader than 
the [C i] line widths which probably reflects changes of the 
optical depths. 

Compared to the northern field, the ^^CO and ^^CO 
lines are weaker while the [C i] 1-0 line is stronger. The 
CO 4-3 peaks at only ~ 20 K. The hue ratios of CO 
7-6/4-3 and [Ci] 2-1/1-0 at the interface position are 
0.3 ± 20% and 0.6 ± 20% respectively. These ratios are 
slightly lower than in the northern field, suggesting lower 
temperatures and densities. However, optical depth effects 
and self-absorption may also play a role. The [C i] 1-0/CO 
4-3 ratio is 0.34 at the interface, a factor 2 higher than 
in the northern interfa ce. The larg e -scale survey of the 
entire Carina region bv lZhang et al.l ( 20011 ) showed varia- 
tions between 0.14 and 0.45, similar to the ratios we find 
at angular resolutions which are a factor ~ 5 higher. Gas 
kinetic temperatures must be at least 30 K to explain the 
^^CO 4-3 main beam temperatures. In the southern in- 
terface, this is consistent with Tox([Ci])=30.5K (^sk) de- 
rived from the [C i] line ratio. 

The LTE analysis at the southern position indicates a 
C/CO abundance ratio of 0.26, slightly higher than the 
abundance ratio in the north. Taking into account the 
calibration errors, the C/CO abundance ratios lie between 
r-^ 0.1 and ~ 0.4 at both positions (Table[T]). Similar abun- 
dance ratios have been observed in other Galactic star 
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forming regions like CepheusB or NGC 7023 (sec refer- 
[Mookeriea et al]l20Q6l) . 



ences m 




500 1000 1500 2000 
Frequency [GHz] 



500 1000 1500 2000 
Frequency [GHz] 



Fig. 9. Integrated intensities at the two interface posi- 
tions in Carina-North and South (Tabled]). Filled symbols 
show the observed ^^CO 2-1, [Ci] 1-0, 2-1, ^^CO 4-3, 
and 7-6 intensities. Error bars denote the 15% ealibra- 
tion error. Model results of the best fitting elumpy PDR 
model (Table[2|) are shown by sohd (^^CO) and dotted 
lines (^'^CO), and open symbols. 



5. Clumpy interface regions 

Figureini shows the cooling intensitieqj observed at the two 
interface positions (Table[T]) along the cuts presented in 
the previous section. At both positions, the 7-6 intensity 
exceeds the 4-3 intensity. As this transition has a criti- 
cal density of 1.4 10^ cm~^ and an upper- level energy cor- 
responding to 157 K, temperatures and densities are ex- 
pected to be similarly high. 

In the following, we compare the observed e mission 
with the predictions of a clumpy PDR model (jCubick 



20051 ICubick et"an l2007t ). We assume that the emis- 



sion stems from an ensemble of spherically symmetric 
clumps within the beam. The clumps follow the canon- 
ical clump mass and mass size distributions: dN/dM ex 
M~" and M ex K^ . For the clump mass spectrum de- 
rived from molecular line observations, a ~ 1.8 in a 



large number of G alactic clouds (e.g. iKramer et al.lll998t 
Simon et al.ll2001 ) . For the power law describing the mass- 
size distribution, 7 ^ 2 .3 in several Galactic clouds (e.g. 
Heithausen et al.l 119981 ). As a result, the clump density 



distribution is n ex M^ '^^'^ , i.e. the density increases with 
decreasing clump mass. For the nearby clouds L 1457 and 



the Polaris Flare presented in the Kramer and Heithausen 
papers, the canonical clump mass and mass size distribu- 
tions continue down to masses of less than 10"'^ Mq and 
radii of less than 10~^ pc. 

The emission from a clump of given m ass and den- 



sity i s calculated by the_ KOSMA-t model (jRolhg et al 
120061 : 

n(r) 



Storzer et al.l Il996l ). The radial density structure 
with radius Rq is given by 



clump 
-1-5 for 0.2 



7'/-Ro < 1 and n{r) 



We used / I^du [erg s"^ cm-2sr-^] = 2|^ / Tdv [K kms^^] 



of a model 
n{r) ex no (r/i?o)"'" for 0.2 < 
11.18 no = const, for r/Ro < 0.2. The clump averaged 
density (n) is about twice the density at the clump sur- 
face, called henceforth surface density ng: (n) ~ 1.91 ng. 
PDR models of individual clumps were calculated on a 
grid of FUV fields log(x) = 0,0.5, to 6, clump surface 
densities log(no/cm^'^) = 2,2.5, to 7, and clump masses 
log(Afci/MQ) = —3,-2.5, to 2. The parameters control- 
ling the ensemble model are the FUV field, the average 
ensemble density, the total ensemble mass, and the lower 
and upper clump mass limits. Here, we assume that all 
clumps are heated by a constant FUV field, as derived 
from the HIRES observations presented above. We chose 
a lower mass cutoff of the ensemble of 10~^ Mq and an up- 
per cutoff of 100 M0 . Below, we show that the exact choice 
of these limits docs however not have a stron g impact on 
the resulting model intensities (cf. ICubick et a l. 2007). To 
fit the observed intensity ratios [Ci] 2-1/1-0 and CO 7- 
6/4-3, the only free parameter is then the average density 
of the clump ensemble (n)ons- 

We adapted the total ensemble mass to fit the absolute 
intensities. The best fitting models and their parameters 
are shown in FigurelH and in Tabled Table [3] lists the ra- 
tios of observed over modelled intensities for all 5 observed 
transitions. For a perfectly fitting model, they would all 
equal 1 ± 15%, refiecting the estimated calibration error 
of the observed data. 



Northern interface. An average ensemble density of 
2 10^ cm~'^ allows us to reproduce the observed ^^CO 7- 
6/4-3 and [Ci] 2-1/1-0 line ratios very well at the north- 
ern interface position. Lower densities do not reproduce 
the high CO 7-6/4-3 line ratios, but rather lead to a peak 
of the cooling curve at J < 5. The absolute intensities of 
the ^^CO and ^^C lines can be reproduced to within 20% 
for an ensemble mass of 400 Mq . The strongest deviation 
occurs for ^^CO; the observed ^'^CO 2-1 line intensity is 
only 40% of the modelled intensity. 

The geometrical beam-filling factor of the clump en- 
semble is defined as the ratio of the sum of all clump 
areas over the beam area, i.e. (pA = ^ens/^boam, with 
f^ens = J2'^^ci/d^ at distance d. The beam-filhng factor 
is 3.9 at the northern interface position (Tabled]). We as- 
sume that clumps don't overlap, i.e. that their emission es- 
capes without being absorbed by foreground clumps, and 
there is no interclump medium. Clumps have a constant 
velocity FWHM of Av = 1.7kms~^. The clu mpy ensem- 
ble models developed by ICubick et al.l ( 2007 ) do not yet 
predict line profiles. Table|2] gives the densities and radii 
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Table 2. Parameters of the clump ensembles at the two positions which fit best the observed intensities. Columns (2) to (6) 
list the input parameters from which the values given in columns (7) to (11) are derived: x is the FUV field in Draine units, 
(n)ens is the average density of the clump ensemble, M™'°, M™*"^, Mtot are the minimum and maximum clump mass, and the 
total ensemble mass, respectively, no.min and no, max are the smallest and largest clump surface densities, 7?min and Rmax are 
the smallest and largest clump radii, and (I)a is the beam-filling factor of the clump ensemble. 



Aa/A(5 
(1) 


X 
(2) 


(n)cns 

cm~^ 
(3) 


Mo 
(4) 


Mo 
(5) 


Mtot 

Mo 

(6) 


no,min 

cm-^ 

(7) 


no, max 

cm^^ 

(8) 


pc 

(9) 


-Km ax 

pc 
(10) 


4>A 

(11) 


-11.7/6.8 
+0.8/ - 7.2 




210" 
210-'^ 


10"^ 
10-2 


10^ 
10^ 


400 
220 


5.310* 
5.310* 


8.8 10" 
8.8 10" 


3.710-^ 
3.710-^ 


2.1 10-" 
0.2 


3.9 
2.1 



Table 3. Ratios of observed integrated intensities over integrated intensities predicted by the clumpy PDR model. The input 
parameters used are given in Table[2l deviations from this setting are listed in column (7). 
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^'^CO 


CO 
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(',') 


2-1 
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7-6 
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(1) 


(2) 


(3) 


(4) 


(5) (6) 


(7) 




-11.7/6.8 


0.4 


1.1 


1.2 


0.8 0.8 






-11.7/6.8 


0.3 


0.9 


0.8 


0.7 0.8 


Mmin 


= 10 -^ Mo 


-11.7/6.8 


0.4 


1.3 


1.7 


0.8 1 


Mmin 


= 10 -'Mo 


-11.7/6.8 


0.5 


1.7 


1.5 


1.4 1.5 


Mmax 


= 10^ Mo 


-11.7/6.8 


0.3 


0.8 


1 


0.5 0.5 


jVJmax 


= 10' Mo 


+0.8/ - 7.2 


0.6 


0.8 


1 


1.2 1.2 







of the smallest and largest clumps of the ensemble. The 
mass of the smallest clump is set to IQ-^Mq. The virial 
mass of the smallest clump, Mvir oc RAv^ is a factor 400 
larger. In other words, the velocity width of these clumps 
are far too large for pure gravitational virialization. The 
external pressure which would be needed to confine these 
clumps, lO^Kcm-^, is far too large to be maintained by 
an interclump medium. The large ratio indicates instead 
that these clumps, if they do exist, arc transient features 
of the turbulent gas and evaporating on very short time 
scales: tovap ~ R/(Tv = i?/FWHM\/8 In 2 w 5000 yrs. 

The modelled intensities depend only weakly on the 
lower and upper mass limits of the ensemble as shown 
in Table[31 Varying one of the mass limits by one order 
of magnitude while keeping all other 4 input parameters 
listed in Table[2] unchanged, the absolute intensities of the 
5 transitions vary by less than a factor of ^^ 2 relative to 
th e best fitting s olutio n. This is discussed in more detail 



Cubick et al. 



( 20071) . 

The run of ^^CO intensities with rotational number J, 
i.e. the CO cooling curve, predicted by the model (Fig. [5]) 
shows a maximum at J = 7 and a steep drop at higher ro- 
tational numbers up to J = 11 by more than a magnitude. 
A secondary peak of the cooling curve is seen at J = 15; it 
is almost a factor 10 weaker than the first maximum. The 
second peak shows up only for FUV fields above ~ 10^ 
and it stems predominantly from the surface regions of 
the small, dense clumps. These clumps are large in num- 
ber, have a large total surface area, which is heated to 
high temperatures, sufficient to excite these high-J lines. 
The abundance of CO is increased in the hot gas of the 
H1/H2 surface layer, near the OH density peak, due to 



secondary producti on paths which become imp ortant at 
high temperatures (jSternberg fc Dalgarndll995l ): 



OH 
OH- 
C0+ 



C+ 

c+ 
-H2 



C0 + 

C0+ 

HCO" 



H 



(R13) 
H^CO + H^^ 

fH (R127) 



(R125,R126) 



HCO+ + e -^ CO + H (R129) 



We suspect that these reactions lead to the secondary peak 
of the cooling curve shown in FigurelHl 

Southern interface. At the southern interface position, the 
FUV field is only 10^'^ xo- An average ensemble density 
of 2 10^ cm"-^ reproduces the observed ^^CO 7-6/4-3 and 
[C i] 2-1/1-0 fine ratios. Within 40%, the clumpy model is 
also consistent with the absolute intensities of all 5 tracers 
for a total ensemble mass of 220 Mq (Fig.El TableslHO]). 
Due to the weaker FUV field compared to the northern 
position, the CO cooling curve peaks at slightly lower fre- 
quencies and a second peak is not visible. Similar to the 
northern position, the strongest deviation occurs for the 
13 CO 2-1 lines. 

6. Summary 

We have mapped the emission of atomic carbon and CO 
in two fields of the GMCs surrounding yyCar. Combining 
CO 4-3 and [Ci] 1-0 data with CO 7-6 and [Ci] 2-1 
data allows us to study the excitation conditions of both 
gas tracers. The northern 4' x 4' field lies adjacent to the 
compact OB cluster Tr 14 and is associated with the H 11 
regions Carl-E and Carl-W. The spectral lines show a 
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rich kincniatical structure, with several velocity compo- 
nents along the lines of sight spread between —10 and 
— SOkms""'^. The southern 4' x 4' field lies on a molecular 
ridge south of rjCsn. Here, spectral lines show only one 
rather narrow Gaussian shaped velocity component. The 
FUV field derived from HIRES/IRAS far-infrared data is 
X ~ 10^'^, about a factor 10 weaker than in the north- 
ern field. IRAC/Spitzer maps at 8/zm show the detailed 
structure of the FUV illuminated dust. The edge of the 
southern region is clearly delineated by the 8^m emission, 
while the northern region shows a more complex spatial 
structure. 

In both regions, CO and [Ci] emission is co-extensive. 
Both species appear to trace the bulk of molecular gas, 
rather than the interface regions. This coincidence has 
been found in many of the previous large-scale studies of 
Galactic clouds c onducted with the Mt.Fuji, KOSMA, and 
other telescopes. iPapadopoulos et al.l (J2004t ) stress the im- 
portance of cosmic rays in raising the C/CO abundance 
ratio and argue that dynamic and non-equilibrium chem- 
istry processes explain the ubiquity of [C i] within molec- 
ular clouds as well as at PDR interfaces. Here, we argue 
that steady state, clumpy PDR models provide an alter- 
native explanation. 

The northern and southern Carina regions both show 
a slight brightening of the CO 7-6 and [C i] 2-1 transitions 
relative to the 4-3 and 1-0 transitions near the interfaces. 
In the south-east part of the northern region, the upper 
transitions brighten towards the interface (Fig.[3|). These 
transitions also show increased intensities near the embed- 
ded IRAS point source in the southern region (Fig. [7]). 

We selected two cuts from the cloud cores through 
the interface regions, roughly pointing towards the illu- 
minating sources. The southern cut lies to the east of 
IRAS 10430-5931, avoiding its complicated structure and 
additional heating sources. The northern cut runs from 
the peak of [C ii] emission (JBrooks et al.l 120031) through 
the H II regions Car I-W and I-E across the edge-on inter- 
face towards the center of the Tr 14 cluster. 

The peak ^^CO line temperatures along the northern 
cut indicate kinetic gas temperatures of at least 50 K, 
while a lower limit of 30 K is found along the southern 
cut. 

We selected two positions at the northern and southern 
interfaces for a detailed analysis. The ^^CO 7-6/4-3 and 
[Ci] 2-1/1-0 intensity ratios at these positions are 0.3-0.4 
and 0.6-0.7, respectively. We use PDR models to interpret 
the observed line intensities. Our models assume that all 
observed emission stems from an ensemble of spherically 
symmetric PDR clumps and there is no emission from 
an interclump medium or a diffuse halo surrounding the 
clouds. PDR clu mps are modelled using the stationary 
KOSMA-r code (|Rollig et al.ll2006h . The clump distribu- 
tions follow the canonical mass and radius distributions 
found for molecular clouds. At the two positions analyzed 
in detail, we find that clumpy PDR models are consis- 
tent with the observed absolute intensities of the ^^CO 



and [Ci] lines to within 20%, i.e. at about the calibration 
accuracy. 

Since the line ratios observed at the two interface posi- 
tions are fairly typical for the entire observed regions, we 
conclude that stationary, clumpy PDR models can simul- 
taneously reproduce the observed [C i] and CO emission 
of the lower transitions and the upper CO 7-6 and [C i] 
2-1 transitions. 

However, the observed ^^CO 2-1 intensities are only 
about half of the modelled intensities at the north- 
ern and southern posit ions, respectively. Interestingly, 
IPineda fc BenschI ( 20071 ) report a similar finding. They 
analyzed "'^^CO and [Ci] emission observed in the dark 
cloud globule B 68. This emission can be reproduced us- 
ing a single KOSMA-r spherically symmetric PDR model, 
while the observed ^^CO intensities are a factor ~ 2 too 
low. 

To study the gas under a broad range of conditions, 
we plan to extend the present [C i] and CO maps to fully 
cover the GMCs surrounding TyCar using the SMART 
492/810 GHz mu ltipixel array receiver at NANTEN2 



(JGraf et al.ll2003[ ). Future observations with APEX and 



Herschel will tell in how far the emission of other species 
tracing the detailed photo-chemical network is also con- 
sistent with stationary, clumpy PDR models. 
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